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O : ABSTRACT 
(N ; 

^ I Ultra-high energy cosmic rays (UHECRs) may originate from the decay of 

^ I massive relic particles in the dark halo of the Galaxy, or they may be produced 

I in nearby galaxies, for example by supermassive black holes in their nuclei. The 

^ I anisotropy in the arrival directions is studied in four dark halo models (cusped, 

if^ I isothermal, triaxial and tilted) and in four galaxy samples (galaxies intrinsically 

^ I brighter than Centaurus A within 50 and 100 Mpc, and galaxies intrinsically 

m ; brighter than M32 within 50 and 100 Mpc). 

2 ' In decaying dark matter models, the amplitude of the anisotropy is controlled 

O ■ by the size of the Galactic halo, while the phase is controlled by the shape. 

' As seen in the northern hemisphere, the amplitude is ~ 0.5 for cusped haloes, 

I ' but falls to ~ 0.3 for isothermal haloes with realistic core radii. The phase 

^ ■ points in the direction of the Galactic Centre, with deviations of ~ 30° possible 

^ ■ for triaxial and tilted haloes. The effect of the halo of M31 is too weak to 

• provide conclusive evidence for the decaying dark matter origin of UHECRs. In 

■ extragalactic models, samples of galaxies brighter than Centaurus A produce 

' substantial anisotropies (~ 1.8), much larger than the limits set by the available 



data. If all galaxies brighter than M32 contribute, then the anisotropy is more 
modest 0.5) and is directed toward mass concentrations in the supergalactic 
plane, like the Virgo cluster. 

Predictions are made for the south station (Malargiie) of the Pierre Auger 
Observatory. If the UHECRs have a Galactic origin, then the phase points to- 
wards the Galactic Centre. If they have an extragalactic origin, then it points in 
the rough direction of the Fornax cluster. This provides a robust discriminant 
between the two theories and requires ~ 350 — 500 events at South Auger. 

Subject headings: Cosmic rays: origin, anisotropy - galactic halo - dark matter 



PACS classification codes: 96.40, 98.70. S, 95.35, 98.35 
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1. Introduction 

Cosmic ray particles with energy in excess of ~ 4 x 10^^ eV have been detected by a 
number of independent experiments over the last few decades . These include the Haverah 
Park 1^ , Fly's Eye and Akeno Giant Air Shower Array (AGASA) |0] experiments, which 
together have recorded over a hundred such ultra-high energy cosmic rays (UHECRs). It 
has long been recognised that the existence of these UHECRs poses an awkward problem [^. 
At such high energies, the typical range of protons decreases rapidly because of interactions 
with the cosmic microwave background photons, becoming as low as ~ 20 Mpc at the highest 
observed energy of ~ 3 x 10^*^ eV It is even smaller for nuclei ||^. There is evidence for 
a change in composition towards protons at the highest energies from the elongation rate of 
the showers observed by Fly's Eye and its successor HIRES as well as AGASA (when 
appropriately reanalysed) . If the UHECRs are nucleons and their sources are distributed 
homogeneously throughout the universe, this should create a "GZK cutoff" in the energy 
spectrum at -Eqzk ~ 4 x 10^^ eV @. This can be evaded if the UHECRs are neutral particles, 
e.g. neutrinos or photons. However, neutrinos which interact only weakly have too small a 
cross-section to initiate the observed airshowers ||10|-Q Photon-initiated showers would have 
different characteristics than those observed; in particular, the highest energy Fly's Eye event 
was probably a nucleon 0, |1^ . A recent analysis of horizontal showers in the Haverah Park 



data sets a bound of 55% on the photon component (and a bound of 30 % on the iron nuclei 
content) of post-GZK UHECRs |]13[. So, if cosmic rays with energies E > Eq2,k <^^^ protons. 



as we assume on the basis of the above evidence 0, ^, |T2[, then they they must originate 
within the Local Supercluster. 

There has been a plethora of suggestions as to the possible origin of UHECRs, which may 
be divided into 'top-down' and 'bottom-up' models. In the former class, the extreme energies 
are provided by the decay of relic topological defects or super-massive particles. In the latter 
class, the extreme energies are provided by the acceleration of particles in astrophysiscal sites 
such as gamma-ray bursts and active galaxies. In this paper, we investigate the observational 
signatures of one 'top-down' and one 'bottom-up' model in some detail. 

One 'top-down' possibility is that the UHECRs may originate from the decay of super- 
massive relic particles which may constitute (a fraction of) the dark matter in galactic haloes 
1^, |T5[. Such particles must have a mass ^ 10^^ GeV to account for the highest energy events. 



with a lifetime exceeding the age of the universe. A well-motivated particle physics candidate 



^If new physics, such as large extra dimensions (implying TeV scale quantum gravity), is invoked, the 
neutrinos may interact more strongly through Kaluza-Klein graviton exchange. It appears unlikely that the 
increased cross-section would be adequate to create the observed air-showers |]lT| 
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with both the required mass and metastabihty, proposed before the definitive Fly's Eye event 
0, is the "crypton" — the analogue of a hadron in the hidden sector of supersymmetry 
breaking in string theories []TB[. It has been recently noted that such particles can be readily 
produced with a cosmologically interesting abundance in the time- varying gravitational field 
at the end of infiation [|l7l . A key test of this model is the expected small anisotropy caused 
by the offset of the Sun from the centre of the Galaxy UTSf. It has been claimed that the 
expected anisotropy is already higher than that observed at lower energies ~ (1 — 5) x 10^^ eV 
|19| . However, this is beside the point, as such measurements do not constrain the new 'flat' 
spectrum component of cosmic rays which extends beyond the GZK cutoff [^]. It has also 
been claimed that the absence of any excess signature in the direction of the halo of the 
nearby Andromeda galaxy is an "insuperable" problem for this model |T9|. This claim has 



already been contested on the grounds that any excess signal would not be detectable in the 
present sample of ~ 0.5 events per 10° x 10° solid angle |2^. A number of authors p!8| - pT | 
have calculated the expected anisotropy assuming either an isothermal halo proflle or the 
"universal" density proflle suggested for cold dark matter dominated halos |22|. However, 



very little is known for certain about the shape, proflle and extent of the dark halo of the 
Galaxy, and so it is important to examine more general halo models. 

One 'bottom-up' possibility is that the UHECRs originate in nearby extragalactic sources. 
Estimates of the intergalactic magnetic fleld (IGMF) are typically of C(10~^) Gauss [P^ , 
which is insufficient to defiect such energetic protons by more than a few degrees. Hence, the 
UHECRs should point back in the direction of the sources. Nevertheless, searches for sources 
within ~ 3° of observed events have not had any success [^. The UHECRs (particularly in 
the Haverah Park sample P]) show a correlation with the supergalactic plane but this is 



not seen in the Fly's Eye p6[ or AGASA data. The latter dataset has the largest number 
of events and a detailed analysis finds the arrival directions are consistent with isotropy, but 
with some evidence of clustering on an angular scale of ~ 2.5° close to the supergalactic 
plane [^. Various possibihties have been suggested for altering the IGMF to enable an 
origin in nearby radio galaxies or active galactic nuclei, e.g., M87 or Centaurus A but 
these have difficulties 



For UHECRs to originate in nearby extragalactic sources, they 
must come from a population that is larger than the nearest active galactic nuclei or radio 
galaxies. It is believed that all big galaxies harbour supermassive (~ 10^ — 10^ M0) black 
holes and have passed through an active phase in the past, even if they are quiescent today 
Processes are known that permit efficient energy extraction from supermassive black 



holes ||3T|, and so one possibility is that nearby galactic nuclei are the arena of production of 



UHECRs even though detailed pathways are lacking. Studies have already shown that 
if the sources of UHECRs trace the density field of nearby galaxies, then the observed fiux 
can possibly be recovered, at least for hard injection spectra p3|. Although the energetics of 



-4- 



the production of the particles may be problematic, this is nonetheless an interesting model 
to contrast with the decaying dark matter hypothesis. It also gives rise to an anisotropy 
signal in the arrival directions of UHECRs, as the nearby rich clusters (like Virgo and For- 
nax) make a disproportionately large contribution. We must keep in mind however that the 
expected anisotropy may be diluted by a factor of up to ~ 2 by the concomitant isotropic 
background from even further sources. In practice, given that this dilution is sensitive to the 
assumed injection spectrum of the sources (e.g., it would be negligible for an spectrum 
I), we do not explicitly take this into account. 



In this paper, we always neglect the effects of magnetic fields on UHECR trajectories. 
This is justifiable if the sources lie in the Galaxy's halo, as studies suggest that the influence 
of the Galactic magnetic field on trajectories is small if the energies are greater than 4 x 
10^^ eV It is a more controversial assumption if the sources are extragalactic. If an 



inhomogenous energy density distribution which follows the Lyman-a forest distribution is 
assumed |^6|], then the estimates of the IGMF can be as large as (9(10^®) Gauss. It also 



possible that the IGMF has a small filling factor and is structured on large scales of up to 
tens of Mpcs |37]. Hence, it is sometimes claimed that a plausible upper limit to the IGMF 



could be as high as O{10~^) Gauss, in which case considerable deflection of incoming charged 
UHECRs is a real possibility . However, there are also strong arguments against such high 
values for the IGMF [^. In particular, it is hard to see how such a field can be generated 
dynamically over tens of Mpcs, as the eddy turnaround time is larger than a Hubble time. 
Note that our choice of 4 x 10^^ eV as the threshold is a reasonable compromise, i.e., it is 
high enough to minimize magnetic deflection without losing too much flux. In principle, 
a lower energy cut may appear attractive for analysing the halo source model, since the 
propagation distance is less than ~ 100 kpc. However, to fit the spectrum still requires 
the new 'flat' component of cosmic rays (from dark matter decays in this case) to remain 
subdominant below the i?"^'® power-law extension from lower energies until ~ 4 x 10^^ eV 
|0[. Moreover, this cut enables easy comparison with earlier studies of the expected 



anisotropy H!], 02 



The motivation for the paper is that the prospects for detection of any anisotropy 
signal in the UHECRs in the next few years are good. The southern station of the Pierre 



Auger Observatory |43| is already under construction at Malargiie, Argentina and will be 
complete by 2004. The array consists of 1600 water-Cerenkov detectors distributed in a 
grid covering about 3000 square kilometers. The northern station is planned for Utah and 
will then enable continuous, full-sky coverage. Estimates have already been made of the 
number of events that Auger will detect after five years of operation [Q. It is reckoned that 
there will be ~ 2200 above 4 x 10^^ eV, ~ 250 above 10^° eV and ~ 35 above the highest 
energy Fly's Eye event (3.2 x 10^° eV). Another experiment that will have a considerable 
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impact is the Extreme Universe Space Observatory EH]. This is scheduled for flight on the 



International Space Station starting in mid 2007. It will detect flourescent light produced 
when UHECRs interact with the Earth's atmosphere. Given the richness of forthcoming 
experimental opportunities, it is interesting to examine possible tests. Can we use the 
anisotropy signal to distinguish between a local origin in the halo of the Galaxy and a more 
distant origin in nearby galaxies? If the UHECRs originate in the Galaxy's halo, what can 
be learnt about its structure? If the UHECRs emanate from nearby astrophysical sources, 
what can be learnt about their distribution? 

The paper is organised as follows. Section 2 discusses the theoretical framework for 
analysis of the anisotropy signal. Section 3 introduces four halo models, while Section 4 
constructs four samples of nearby galaxies. The prospects for distinguishing between the 
galactic and extragalactic origins of UHECRs with the existing Haverah Park and AGASA 
experiments and especially the southern station of the Auger Observatory are examined. 
Section 5 sums up our conclusions. 



2. The Anisotropy Signal 

2.1. The Detector Response 

The emissivity of UHECR per unit volume is proportional to the number density of 
sources n{r). The incoming flux of UHECR per unit solid angle F as a function of right 
ascension a and declination 6 is 

F{a,6)(x / (isn(R0 + s(q;, 5)), (1) 

Jo 

where s is the heliocentric position vector. Here, Smax marks the extent of the distribution 
of sources, which may depend on direction. 

The measured flux of UHECR is the incoming flux modulated by the response of the 
detector h{6), which is the relative efficiency of the detection of events with direction. It 
depends only on the declination and not on the right ascension, because, if the detector is 
run with reasonable efficiency, there will be almost uniform exposure in right ascension after 
a year or so. The response as a function of declination depends partly on the attenuation of 
showers in the atmosphere, partly on the sensitivity of the detector as a function of zenith 
angle. 

For the scintillator arrays of AGASA and the water- Cerenkov detectors used at Haverah 
Park, the declination distribution h{S) of the observed events is reproduced in Figure 1 of 
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Uchihori et al. Our procedure is to fit third-degree polynomials (conveniently bounded 
by step functions) to these experimental curves. For AGASA, this gives 

hAciS) = 0.323616 + 0.0361515 5 - 5.04019 x lO"'* 5^ + 5.53941 x 10"^ 5^ (2) 

valid for —8.0° < 6 < 87.5°. For Haverah Park, this gives 

h^p{5) = 0.108234 + 0.0320955 5 - 1.881358 x 10"^ 5^ - 1.948817 x 10"^ 5^ (3) 

valid for 6 > —3.3°. These expressions are normalised to a maximum value of unity. 

In the absence of measurements from South Auger, we explore two possibilities. First, 



following Medina Tanco & Watson we use the declination distribution found at Haverah 
Park (54° N), mirrored and shifted to the latitude of the South Auger site at Malargiie, 
Argentina (35° S). The motivation for this is that the peak of the response function usually 
coincides with the geographic latitude (at least when there is not continuous exposure at the 
Poles). So, we use: 

hi5) = hnp{19°-6). (4) 
Second, following Sommers ^M, we use the analytic expression valid for a detector at a 



single site with continuous operation and constant exposure in right ascension. Suppose the 
detector is at latitude A and that it is fully efficient for particles arriving with zenith angles 
6 less than some maximum value 6'm, which is taken as 60°. This gives us: 

/i2(5) oc cos A cos 5 sin + ttm sin A sin (5, (5) 

where is the local hour angle at which the zenith angle becomes equal to 6*111- It is given 
by 

[0 ife>l 

ttm = < TT if ^ < -1 (6) 

[ acos(,^) otherwise 

and 

cos 6'iii — sin A sin S 
COS A cos 

Note that the two declination distributions hi{6) and h2{S) are quite different. The maximum 
of h2{S) is (generally) in the direction of the Pole, whereas the maximum of hi{6) is in the 
direction of the geographic latitude of the South Auger site. Often, but not always, our 
results do not depend too much on the details of the response function. 
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Fig. 1. — Plots of the prob bility distributions -P(0; where ^ = 5/7^, (left panel) and P{{}) 
(right panel) as given by eqns (|12D and ([T3|) respectively. The curves are labelled with the 
relevant value of fco = \7i?N . Directly measurable quantities are only involved through this 
quantity. 



Experiment 


Amplitude 
7^ 


Phase 


N 


4 

(95%) 


^5 

(95%) 


(95%) 


(95%) 


AGASA 


0.262 


298.5° 


57 


oo 


0.209 


0.242 


145° 


Haverah Park 


0.583 


173.6° 


27 


0.622 


0.452 


0.491 


95° 


Auger I 


0.262 


298.5° 


1000 


0.09 


0.09 


0.09 


19.6° 


Auger II 


0.583 


173.6° 


1000 


0.09 


0.09 


0.09 


6.5° 



Table 1: Summary of the experimental data for AGASA and Haverah Park, together with 
two projections for Auger. The first projection assumes that the amplitude and phase are 
the same as that seen by AGASA, the second the same as that seen by Haverah Park; in 
both cases the number of events is increased to 1000. The amplitude TZ and phase ip are 
calculated from the number of UHECRs, together with the 95% confidence limits. (For 
AGASA, the 95% confidence limit is not strictly speaking well-defined, as the integration 
from 7?. to oo gives a probability of only 0.449) 
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0% 20% 40% 60% 80% 100% 



Fig. 2. — Contour plots of the UHECR flux in equatorial coordinates for our four dark halo 
models, namely (a) cusped, (b) isothermal, (c) triaxial and (d) tilted. These are Hamer- 
Aitoff projections. The Galactic plane is marked in each figure. The effect of the halo of 
M31 is visible in the upper right of each plot. 




Fig. 3. — As Figure H, but now the UHECR flux has been convolved with the response in 
declination for the southern site of the Auger Observatory. 
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2.2. Harmonic Analysis 



The strength and direction of the anisotropy in the arrival directions of the UHECRs 
can be quantified by harmonic analysis, as originally suggested by Linsley |^ . The number 
of events is proportional to N, where: 



N = J d6 J da cosSh{5)F{a,S). 



Let us define 



2 

" = N 

so that the amplitude TZ and phase ip of the first harmonic is 



d6 J da cos 6 cos a h{6)F {a, 6), 
d6 / da cos 6 sin a h{6)F {a, 5), 



ip = atan ( ~ j • 



(9) 



(10) 



Once we have a set of experimental UHECRs with amplitude TZ and phase iJj, the probability 
P{S, i)) that this data came from a population characterized with an amplitude in {S,S+dS) 
and phase {i), i) + di)) is given by 



N 



27r3/2/o(A;o/2; 



exp 



N 



(52 + i7^2-27^5cos(^9-V) 



'11^ 



where = ^N7i? and Jq denotes the zeroth order Bessel function of the first kind. Defining 
^ = S/TZ, then the differential probability distributions P{C) and P{d) are obtained by 
straightforward integration as 



p{^) 



2t 



1/2 



exp 



^(l + 2a]/o(2fcoe), 



1 ±erf 



i/^y^ cos('i9 — ip) 



exp 



/co(cos^('(9 — -ip) 



(12) 
(13) 



7rV2/o(A:o/2) 
1 

27r/o(fco/2) 

Figure |I| shows these probability distributions. For large numbers of events, these distribu- 
tions tend towards normal. It is not obvious that the assumption of normal distributions is 
appropriate for the AGASA and Haverah Park datasets, as the number of events is rather 
modest. In this paper, we always calculate the dispersions about the expected value corre- 
sponding to a given confidence limit by solving the integral equation over the distribution. 
For example, the dispersion corresponding to the 95% confidence limit is given by 



P{^)d^ = 0.95. 



(14) 
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The symmetric dispersion on the amphtude ag can be defined in the same way. However, the 
confidence hmits on the amphtude S are asymmetric about the average values. So, sometimes 
it is useful to use asymmetric dispersions and corresponding to 95% confidence limits 
defined via 

P(nd^ = 0.475, (15) 

(€> 
and 

P{Od^ = 0.475. (16) 

Unless stated to the contrary, we use the asymmetric definitions of confidence limits for the 
amplitude. This distinction can be important at small /cq- 



2.3. The Data 

Table |I| summarises the available experimental data. The Haverah Park shower array 
detected 27 UHECRs with energies above 4 x 10^^ eV and zenith angle less than 45° |[47|| . 
For the AG AS A experiment, we use the data on the 57 UHECRs detected before May 2000 
27| . Table records the amplitude TZ and phase as determined from the two datasetsQ. 



The 95% confidence limits on the amplitude and phase are given, from which it is already 
evident that the presently available data are too sparse to claim a secure detection of any 
anisotropy. This is in accord with AGASA's own analysis of the 47 events they observed 
before August 1998, in which they assert that there is no evidence for significant large-scale 
anisotropy on the celestial sphere. 

Sommers has estimated that the likely size of the dataset after 5 years of operation 
of the Auger Observatory is ~ 2200 events. This though is the combined total from both 
North and South Auger Observatories. In this paper, we are interested in the prospects over 
the next few years and so we restrict our analysis to South Auger. To give a feel for what 
it may achieve should there be a real anisotropy in the UHECR sky, we present two further 
examples in Table In the first (henceforth Auger I), we assume that its amplitude and 
phase are the same as that indicated by AGASA, but that the number of events detected 
over ~ 4 X 10^^ eV is 1000. In the second (Auger II), we assume the amplitude and phase 
are the same as that indicated by Haverah Park (which uses water-Cerenkov detectors of 
the same depth as Auger) and the number of post-GZK events is 1000. In both cases, the 



^ Since galactic magnetic fields may influence the trajectories of cosmic rays with energies < 4 x 10^^ eV 
55[, we do not consider these events in the analysis. 
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size of even the 95% confidence limits is reduced to typically ~ 10° in phase and ~ 30% in 
amplitude. 

Table || gives the details of further tests for isotropy. We list the probability that an 
amplitude higher than that measured TZ can arise from an underlying isotropic distribution. 



This is given by the formula [it 



P(>7^) =exp(-A:o). (17) 

We also give the probability that the data comes from a model with an amplitude S between 
(isotropic) and the symmetric dispersion as corresponding to the 68% and 95% confi- 
dence limits. This table underlines that the present AGASA and Haverah Park datasets 
are inadequate for the unambiguous identification of any anisotropy signal, if it exists. For 
example, an amplitude equal to or exceeding that reported by AGASA can come from an 
isotropic distribution with a probability exceeding a third, which is hardly negligible. The 
final two rows of Table ^ show the impact that South Auger will make. The small numbers 
for the probabilities demonstrate that genuine anisotropy signals with the same amplitude 
and phase as indicated by AGASA and Haverah Park will be convincingly detected. 



3. The Dark Halo of the Galaxy 
3.1. Models 

There is little definite evidence on the structure of the Galaxy's dark matter halo. The 
HI gas rotation curve of the Milky Way cannot be traced much beyond ~ 20 kpc. It is the 
kinematics of the distant satellite galaxies that provide the best evidence on the total mass 
and the extent of the Galaxy. The most recent estimate finds a total mass ~ 2 x lO^^M© 



and an extent of ;^ 200 kpc. Such a large extent is in accord with earlier studies of the Milky 



Way and external spiral galaxies [49 



N-body computer simulations of structure formation in hierarchical merging cosmogo- 
nies have provided evidence for a cusped density profile for dark haloes 



n r oc 



(r + re)(r + r^y 



where r is the position with respect to the Galactic Centre and r = |r| is the spherical 
polar radius. Here, Vs is the scale radius, which is typically ~ 10 kpc for the Galaxy, while 
is set by the resolution limit of the simulations, nowadays ~ 0.5 kpc. These models 
(hereafter NFW) have been much in vogue over the last few years. For example, they have 
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already been used in calculations of the flux of UHECRs and the strength of the anisotropy 
signal [0, |2^. The cusped density distribution of the NFW proflle necessarily implies that 



the inner regions are dominated by dark matter and gives a substantial anisotropy signal. 
However, there are three reasons why the NFW is certainly incorrect for the Galaxy ||50|| ; 
flrst, the mass density implied by the luminous disk and bar is already sufficient to account 



for the rotation curve in the inner Galaxy without any contribution from dark matter [51 



second, numerical simulations of barred galaxies shows that copious amounts of dark matter 



in the inner parts slow down and rapidly dissolve bars through dynamical friction and 
third, the microlensing optical depth to the red clump stars already shows that almost all 
the density in the inner Galaxy must be in the form of compact objects capable of causing 



microlensing and so cannot be particle dark matter [B3 



It is more likely that the Galaxy's dark halo has an isothermal distribution with a large 
core radius. Unlike the NFW proflle, this at least satisfles one of the few pieces of secure 
evidence — namely that luminous matter dominates the central regions, whereas dark matter 
dominates the outer parts. As a prototype of this model, we take the cored isothermal-like 
sphere 

(r^ _|_ y-^^Y 

Here, r^. is the core radius, which is chosen to be ~ 10 kpc. This parameter is not constrained 
at all by the observations, but its precise value does not have a big effect on the anisotropy 
of the UHECRs. The model is the spherical limit of a more general family of triaxial density 



distributions [55 



Ax' + By' + Cz' + D 
where r = (x, y, z) and the constants A, B, C and D are 

A = (p-2+g-2_i)^ B = p-\l-p-^+q-^), C = q-\l+p^^-q-^), D = rl{l+p-^+q-^). 

(21) 

Here, p and q are axis ratios of the potential. When p = 1, the halo is oblate, whilst when 
p = q, the halo is prolate. The shape of the Galaxy's dark halo is uncertain. The flaring of 
the outer HI gas disk does provide direct evidence on the vertical force at the midplane from 
the halo, which is greater for flattened models. Analyses of this effect have been used to 



determine the shape of the halo [^, but they are hampered by the unknown contribution 
to pressure support from magnetic flelds and cosmic rays at large distances from the Sun. 
There is also evidence on the shape of the halo from stellar kinematics [^, although this 



is beset with modelling uncertainties as to the orientation of the velocity ellipsoid. Such 
evidence as there is suggests a flattened dark halo with an ellipticity less than E7. However, 
even though evidence on the shape of the dark halo is sparse, it is natural to expect halos 
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built from merging and accretion to possess flattened, generically triaxial, shapes rather than 
spherical ones |0. For the sake of definiteness, we take p = 0.9 and q = 0.75, which means 



that the density contours have semi-axes in the ratio 1 : 0.788 : 0.428. This is a highly 
flattened model with an ellipticity Q of roughly E6. At first glance, a triaxial halo may seem 
to afford a great deal of freedom as to the location of the Sun. This is of course not the 
case! The Sun must be located close to one of the principal otherwise there would 

be gross differences between the Galaxy as viewed at positive and negative longitudes. Both 
local and global dynamical evidence favour the location of the Sun as close to (within 10° 



of) the major axis of the Galactic disk |[59 



There is one other effect that may be common. It has been emphasised |^D[ that the 
outer parts of the dark halo may be misaligned with the disk. This is one possible cause of 
the phenomenon of warping of the neutral gas disk, an effect that is known to be present in 
the outer Milky Way. A family of tilted halo models has density 

where R'^ = x'^ + y'^. The coordinates {x', y', z') are related to (x, z) by a rotation through 
an angle 9 about the x-axis, on which the Sun lies. The motivation for this is that the Sun 
lies nearly on the line of nodes of the warp. We take 9 = 30°, an extreme value, as we are 
interested in the qualitative effects of tilted haloes. 

There are two astrophysical parameters we will need to calculate signatures of the UHE- 
CRs, assuming they arise from the decay of heavy relics in the halo. The first is R0, the 
distance from the Sun to the Galactic Centre, which is taken as (8.5,0,0) The second 
is Tniax; the extcut of the halo, which is taken as ~ 250 kpc in accord with the most recent 
estimates ||8, i9]. The halo extent is poorly known and results depending sensitively on this 



parameter must be viewed with some caution. 



3.2. Results 

Figure ^ shows what our four halo models look like in the cosmic ray sky. These are 
contours of the incoming UHECR flux in equatorial coordinates. Figure ^ shows the contours 
convolved with the response function hi{5) for South Auger. (The figure for h2{6) is very 
similar). If the UHECRs do originate in dark haloes, then the halo of the nearby Andromeda 



•^Recall that the clhpticity is 10 x (1 — 5/a), where b and a are the projected minor and major axes 
respectively. A spherical galaxy is EO, while the most flattened elliptical galaxies are E7. 
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Fig. 4. — Plots of the amplitude against the phase for the four models (cusped, isothermal, 
triaxial and tilted). Also shown are the 95% confidence limits for (a) the AGASA and (b) 
the Haverah Park experiments. The lower panels show the expected impact of South Auger. 
Marked are the 95% confidence limits assuming the amplitude and phase of (c) the AGASA 
and (d) the Haverah park experiments, but with 1000 events. 
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galaxy (M31) may give a tell-tale signature. M31 is the only other large galaxy in the Local 
Group, and it is roughly as massive as the Galaxy [^]. The halo of M31 is modelled as an 
isothermal sphere centered 770 kpc away and with an extent of 250 kpc. In Figures ^ and 
|, the hotspot at ao ^ OOHS'^, 6o ^ 41° (referred to the J2000.0 epoch) is caused by M31. 
It has already been claimed that the absence of a hotspot in the data rules out the decaying 



dark matter origin of UHECRs |jT9[ , though this has been contested by others 



The issue is worth considering in more detail. The ratio of the mass of M31 to that 



of the Galaxy used previously |jT9|, is two to one. There is little dynamical evidence for 
such a ratio. The asymptotic rotation curve of M31 is just ~ 10% higher than that of the 
Galaxy. The most reliable estimate for the mass of the M31 halo based on the motions of 
satellite galaxies, planetary nebulae and globular clusters is ~ 1 x 10^^ M©, though there is 
at least a factor of two uncertainty in this value due to the small sizes of the tracer datasets 
||6^ . On balance, it is probably more accurate to assume that M31 is roughly as massive as 
the Galaxy. Now, to compute the ratio of the total flux received from M31 to that received 



from the Milky Way $m3i/'^'mw, the following expression has sometimes been used pO 



$M3i e Iv^,,nM3idV ^ ^^^^ 



<I>MW D"^ Ivsn'^Mw/r'^ dV 

Here, ^ is the ratio between the masses of the two haloes, D ~ 770 kpc, Vmsi is the volume of 
the halo of M31 and Vsq, is the volume defined by the cone of solid angle 6VL pointing towards 
M31. There are two drawbacks to this formula. First, it assumes that M31 is point-like so 
that nM3i/r^ is taken as constant and equal to its value at the center of M31. Second, it 
assumes that the entire halo of M31 fits into the angular window. Accordingly, we prefer to 
use the exact expression 

$M3i Q-ta/l da fs'-M/i cos 5 dS /;_+ nM3i ds 
^ r;+tt da cos 6 d6 nMw ds' 

Here, the flux ratio is measured in an angular window Aa x A6 centered on the right ascension 
ao and declination 6o of M31. The quantities s± are computed by finding the intersection of 
the line of sight towards with a sphere of radius 250 kpc, while Smax is the usual limit of the 
Galaxy halo. Both the Galaxy and M31 are modelled with identical cored isothermal haloes. 
Table ^ gives the results for the flux ratio computed using the approximate formula (|2^) and 
the exact one (|^). Results are presented when both the Galaxy and M31 are modelled with 
isothermal and cusped haloes of extents of 250 kpc and 100 kpc. We see that the effects 
of an extended dark halo are not always well-represented by the point-like approximation. 
However, formula ( p3D is still valuable, as it gives an upper limit to M31's contribution. 
Given that the expected spread caused by the deflection of a 4 x 10^^ eV proton originating 
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from Andromeda is of the order of a few degrees, a field of view of 10° x 10° is probably 
the most appropriate one to look for the enhancement effect of the M31 halo ||20|. On such 
scales, no substantial enhancement effect is expected and the dominant source of UHECRs 
remains the Galaxy. We therefore agree with the assessment of Medina Tanco & Watson |^ 
- and disagree with Benson et al. |l^ - as regards the visibility of M31 in the UHECR sky. 



Table ^ examines the four halo models (cusped, isothermal, triaxial and tilted). For the 
moment, we assume the canonical values of the halo parameters given in Section 3.1. For each 
model, we compute the theoretical amplitude S and phase "i? in the first two columns. Then, 
the next columns give the probability that the experimental data came from a particular 
model with S in {S — as,S + a^) and in {-d — a^,{} + a^). This is found by computing the 
integral 

P{as,cr^)= _ / PiS,^)dSd^. (25) 

The final two columns give the equivalent results for the one-dimensional distributions of 
amplitude and phase. 

For both AGASA and Haverah Park response functions, the amplitude is greatest for 
the cusped NFW model. This makes good sense, as this model has the largest central 
density. We have argued above that the Galaxy almost certainly does not have such a cusped 
density distribution for its dark halo, and so a better guide to the expected amphtude is 
provided by the other three halo models. Isothermal- like models with a core radius ~ 10 
kpc have an anisotropy of amplitude S ^ 0.4 as measured by the AGASA and Haverah Park 
experiments. For all four models, the phase points roughly towards the Galactic Centre 
at a; ~ 17*^45^,5 ~ —29°. The largest deviation from this direction of 35° occurs for the 
triaxial halo. Notice that the probabilities listed in the final three columns are all far from 
insignificant. The data are compatible with any of the models. The present dataset is 
insufficient to discriminate between possible dark matter halo densities. 

In Table ^ (a), which corresponds to Auger I, we see the improvement in discrimination. 
The amplitude and phase is the same as that seen by AGASA (see Table 0) and happens 
to lie close to the prediction of the triaxial model. Hence, Table § (a) shows that all the 
probabilities are very small, except that of the triaxial model. In other words, there is now 
successful discrimination between halo models. All the probabilities in Table |^ (b), which 
corresponds to Auger II, are very small. The phase, assumed to be that of Haverah Park, 
does not lie close to any of the halo models which are accordingly ruled out. In Table ^ (c), 
we give the expected amplitude and phase for our four halo models and convolved with the 
two different response functions hi{5) and h2{5) proposed for South Auger. These results are 
shown pictorially in Figure ^ The panels show the datasets from AGASA, Haverah Park, 
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Auger I and II, together with 95% confidence hmits. 

The anisotropy of all four halo models has a larger amplitude {S ^ 0.5) at the South 
Auger detector (35° S) than at AGASA (36° N) or Haverah Park (54° N). For the spherical 
models, the phase marks the direction of the Galactic Centre. A signature of triaxiality is 
that the phase has different angular offsets from the Galactic Centre direction at different 
detector locations. 

Figure | investigates the amplitude S and phase seen by South Auger as a function 
of the model parameters for the four halo models. Panels (a) and (b) show the effects of 
varying the scale radius Vg of the cusped NFW halo and the core radius Vc of the isothermal 
halo. These parameters affect only the amplitude of the harmonic and not the phase, which is 
always in the direction of the Galactic Centre. Note that there is no simple way of measuring 
the core radius of a dark halo from astrophysical data. Usually, the core radius is inferred 
by least-square fitting to the galaxy rotation curve, but this is an uncertain procedure as a 
single rotation curve must be decomposed into separate contributions from the bulge, disk 
and halo. If UHECRs are indeed produced by decaying dark matter, these plots raise the 
exciting possibility that halo core radius may become directly accessible to measurement. 
During the first 5 years of operation of South Auger, the core radius could be measured if 
Tc <^ 10 kpc. Panel (c) shows the effect of varying the shape of the triaxial halo. The oblate 
models {p = 1) show largest deviation in phase (~ 15°) from the spherical model {p = q = 1). 
By comparison, stretching or flattening the halo does not change the value of the amplitude 
S by much. In panel (d), the tilt angle is varied. This can give large deviations in phase 
(~ 25°) and significant changes in amplitude. Overall, Figure | suggests the shape of the 
halo controls the phase whereas the scale of the halo controls the amplitude S. 

Figure ^ shows the the number of events A/5 required for a convincing 5a detection of 
the anisotropy signal for the four halo models at South Auger. An approximate method is 
to use 

ATs - (26) 

This gives the number of events required for a signal-to-noise ratio of n„ standard deviations 
in amplitude. This assumes that the probability density P{S) is a normal distribution about 
S = TZ with dispersion equal to \J2/Ms- This is not the case, especially when the number 
of events is small or the amplitude is small (i.e., is small). In an exact calculation, the 
signal-to-noise is computed from its definition 



(27) 
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where 



If we fix n a = 5, as is required for a convincing detection, then eq. ( pTf ) is an imphcit 
equation for Afs, which can be solved iteratively. Figure ^ (a) and (b) show the number of 
events for a 5a detection calculated using the exact expression (full curves) and the normal 
approximation (dashed curves), from which it is evident that the approximation tends to 
underestimate the required number of events. Nonetheless, given the ease with which the 
normal approximation can be computed, it is reassuring to see that it is never seriously in 
error. Roughly 40 events may be needed if the Galaxy's halo has a cusped NFW profile with 
a small scale length rg < 2 kpc. More likely, if the halo is an isothermal sphere with large 
core radii (rc ~ 10 — 20 kpc), then typically 150-500 events may be required. Thus, within 
5 years of operation. South Auger will be able to rule some halo models out. 

Panels (c) and (d) show that shape parameters like triaxiality and tilt do not have a 
substantial influence on the number of required events for a convincing measurement of the 
amplitude. For these models, it is more interesting to ask how many events N^p are needed 
for the error bars on the phase to be reduced to ±/3. Using the normal approximation, this 
gives 

K - (28) 
sm p 

For example, for the 3a error bars on the phase measured at South Auger to be reduced to 
±10° requires ~ 2000 events. This is the number needed to begin to distinguish between 
different triaxial or tilted models and is possible after ^ 10 years of operation. A less 
ambitious target is to ask for the 3a error bars to be ±20°, which requires ~ 500 events. 
This is ample to verify that the phase points in the rough direction of the Galactic Centre 
and thus implicate the decaying dark matter models as the correct mechanism for UHECR 
production. 



4. Nearby Extragalactic Sources 
4.1. Models 

There are two serious impediments to the proposal that extragalactic sources supply the 
UHECRs. First, most of the UHECRs do not point back towards any obvious astrophysical 
sources, such as the active or interacting galaxies within ~ 50 Mpc. For example, the 
possible source of the highest energy Fly's Eye event was examined in and the best 
candidate amongst nearby sources was the starburst galaxy M82, some 37° away. Second, 
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Fig. 5. — Amplitude-phase plots showing how the anisotropy signal measured by South 
Auger varies with the parameters in the (a) cusped, (b) isothermal, (c) triaxial and (d) 
tilted halo models. The legend in (c) and (d) corresponds to the flattening in the potential 
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Fig. 6. — Number of events required for 5a detection of an anisotropy by South Auger for 
the (a) cusped, (b) isothermal, (c) triaxial and (d) tihed halo models. The full and dashed 
curves in (a) and (b) refer to the exact result and the normal approximation respectively. 
The legend in (c) and (d) corresponds to the flattening in the potential q. 
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Experiment 


P(>7^) 


P{Q<S<as) 
(68%) 


P{<d<S<(ys) 
(95%) 


AGASA 


0.377 


0.403 


0.578 


Haverah Park 


0.101 


0.242 


0.494 


Auger I & II 


< 10-3 


< 10-3 


< 10-3 



Table 2: Tests for isotropy on the AGASA and Haverah Park datasets, as well as the pro- 
jected Auger datasets. The first column gives the probability that an isotropic distribution 
could give an amplitude greater than that observed TZ. The second two columns give the 
probability that the data comes from a model with theoretical amplitude S between isotropy 
and the 68% or 95% symmetric confidence interval. For Auger I & II, the probabilities are 
< 10~3^ enabling isotropy to be securely ruled out. 




(d) 



Fig. 7. — Hammer- Aitoff projections in equatorial coordinates-s of the four samples of galax- 
ies, namely (a) Sample I, (b) Sample II, (c) Sample HI and (d) Sample IV. These plots 
record the positions of the galaxies, not the fiux of incoming UHECRs. The Galactic Center 
is marked with a star. 
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(a) Isothermal (Extent of 250 kpc) 



Field Size 


'^'msi/'^'mw 
(Approx, eq. 


23 


) 


'^'msi/^mw 
(Exact, eq. 24) 


2° X 2° 


39.7 


2.5 


10° X 10° 


1.6 


0.6 


15° X 15° 


0.7 


0.4 


(b) Isothermal (Extent of 100 kpc) 


Field Size 


'^'msi/'^'mw 
(Approx, eq. 


23 


) 


'^'msi/'^'mw 
(Exact, eq. H) 


2° X 2° 


16.4 


2.6 


10° X 10° 


0.7 


0.5 


15° X 15° 


0.3 


0.3 


(c) Cusped (Extent of 250 kpc) 


Field Size 


'^*M3i/'^'mW 

(Approx, eq. 


'^'msi/^mw 
(Exact, eq. |§ 


2° X 2° 


18.1 


3.0 


10° X 10° 


0.7 


0.4 


15° X 15° 


0.3 


0.2 


(d) Cusped (Extent of 100 kpc) 


Field Size 


'^'msi/'^'mw 
(Approx, eq. ^3|) 


(Exact, eq. ^ 


2° X 2° 


11.8 


3.0 


10° X 10° 


0.5 


0.4 


15° X 15° 


0.2 


0.2 



Table 3: Ratio between incoming UHECR flux from the M31 halo to that from the Galaxy 
for fields of view of various sizes centered on M31. The second column shows the dangers of 
the point mass approximation, which overestimates the importance of the effect, while the 
third column gives the exact result. In tables (a) and (b), both galaxies are modelled with 
isothermal haloes, but with an extent 250 kpc and 100 kpc respectively. Tables (c) and (d) 
are the same but for the use of the cusped halo models, rather than isothermals. 
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(a) AGASA 



Model 


Amplitude 


Phase 






PM 




S 




(95%) 


(95%) 


(95%) 


Cusped 


0.497 


265.8 


0.328 


0.330 


0.944 


Isothermal 


0.385 


265.8 


0.563 


0.572 


0.944 


Triaxial 


0.251 


295.0 


0.865 


0.898 


0.950 


Tilted 


0.346 


267.8 


0.651 


0.665 


0.945 


(b) Haverah Pcirk 


Model 


Amplitude 


Phase 




Pi^s) 


PM 




S 


1? 


(95%) 


(95%) 


(95%) 


Cusped 


0.446 


265.8 


0.530 


0.977 


0.543 


Isotlioniial 


0.350 


265.8 


0.515 


0.951 


0.512 


Triaxial 


0.224 


301.2 


0.172 


0.887 


0.180 


Tilted 


0.314 


268.2 


0.480 


0.937 


0.513 



Table 4: For each of the models, the theoretical values of the amphtude and the phase are 
calculated using the response functions for (a) AGASA and (b) Haverah Park. We also 
give the probabilities that the experimental amplitude and phase come from within the 95% 
confidence limits of the model. 
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(a) Auger I 



Model 






PM 




(95%) 


(95%) 


(95%) 


Cusped 


< 10-=^ 


< 10-3 


0.095 


Isothermal 


0.014 


0.202 


0.095 


Triaxial 


0.886 


0.950 


0.932 


Tilted 


0.059 


0.509 


0.133 


(b) Auger II 


Model 






PM 




(95%) 


(95%) 


(95%) 


Cusped 


< 10-3 


0.156 


< 10-3 


Isothermal 


< 10-3 


0.001 


< 10-3 



(c) Predictions for South Auger 



Model 


Amplitude 


Phase 


Amplitude 


Phase 




ih) 


ih) 


ih) 


{h2) 


Cusped 


0.883 


265.8 


0.806 


265.8 


Isothermal 


0.558 


265.8 


0.508 


265.8 


Triaxial 


0.551 


257.9 


0.508 


255.9 


Tilted 


0.512 


265.0 


0.457 


264.6 



Table 5: For each of the models, this gives an idea of what we will learn from the southern 
station of the Auger Observatory. Table (a) assumes the amplitude and phase of the AGASA 
experiment, table (b) the amphtude and phase from Haverah Park; in both cases, 1000 events 
are analysed as typical of the number of UHECRs recorded by South Auger. In Table (b), 
only the results for the isothermal and cusped models are recorded; for the other models, all 
the probabilities are less than 10-3. Table (c) presents the predictions for our halo models 
using the two possible response functions hi{S) and h2{S). 



- 26 - 



the possibilities for acceleration of cosmic rays to such extreme energies in any extragalactic 
object seem to be intrinsically limited |6^. Despite the obstacles, the idea that UHECRs 
originate in extragalactic sources is still worth examining. At energies exceeding ~ 4 x 10^^ 
eV, the UHECRs in the Haverah Park dataset exhibit a correlation with the direction of the 



supergalactic plane though this is not seen in the AGASA experiment [^]. However, 
two of the AGASA clusters (CI and C2) do lie close to the supergalactic plane and do have 
tentative identifications with nearby galaxies. 

Let us suppose that supermassive black holes in galactic nuclei produce UHECRs. There 
is increasing evidence that many, perhaps all, galaxies contain central black holes, even if 
their nuclei display only low levels of activity today. Very convincing candidates are provided 
by the nearby galaxy NGC 4258, where maser emission has traced out the gas cloud motions 
around a central black hole of ~ 3.6 x 10^ M©, as well as the Galaxy, where the proper 
motions of stars very close to the center strongly suggest a black hole of mass ~ 2 x 10^ Mq 
Stellar kinematics can now be obtained at high spatial resolution (down to 1") with 



the new generation of spectrographs, enabling the nuclear regions of nearby galaxies to 
be probed. This has strengthened the evidence for supermassive black holes in a number of 
nearby galaxies, such as M31, M32 and M87 The case of M32 is particularly interesting. 



Although it is only a low luminosity galaxy, there is strong evidence for a supermassive black 
hole of ~ 2 X 10® M0 Only its relative closeness enables the black hole to be identified. 



The Third Reference Catalogue of Bright Galaxies (RC3) is reasonably complete for 
nearby galaxies ||6^. It aims to include all galaxies with total B-band magnitudes i?T brighter 
than ~ 15.5 and with a redshift less than 15000 kms~^. We use it to construct the following 
datasets. Samples I and // contain all galaxies in the RC3 intrinsically brighter than M32 
(NGC 221) and closer than 50 Mpc and 100 Mpc respectively. M32 is the smallest galaxy 
for which there is good evidence for a central black hole, so it is a reasonable supposition 
that all galaxies larger than M32 contain black holes, even if quiescent. Samples ///and IV 
contains all galaxies in the RC3 intrinsically brighter than Centaurus A (NGC 5128) and 
closer than 50 Mpc and 100 Mpc respectively. Centaurus A is the nearest active galaxy, so 
these Samples contain only the brightest and largest of the nearby galaxies. 

The motivation behind these choices is to answer the following questions. What is the 
expected anisotropy signal if only big, bright galaxies supply the UHECRs, or if all nearby 
galaxies contribute equally? Is it only galaxies within 50 Mpc that dominate the signal? 
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(a) AGASA 



Model 


Amplitude 


Phase 






PM 




S 




(95%) 


(95%) 


(95%) 


Sample I 


0.309 


111.1° 


0.343 


0.755 


0.523 


Sample II 


0.307 


110.9° 


0.346 


0.759 


0.524 


Sample III 


1.902 


10.7° 


< 10-3 


< 10-3 


0.910 


Sample IV 


1.862 


10.9° 


< 10-3 


< 10-3 


0.910 


(b) Haverah Pcirk 


Model 


Amplitude 


Phase 




PM 


PM 




S 




(95%) 


(95%) 


(95%) 


Sample I 


0.353 


110.8° 


0.783 


0.952 


0.828 


Saiiii)lo II 


0.352 


110.7° 


0.782 


0.952 


0.828 


Sample III 


1.922 


10.8° 


< 10-3 


< 10-3 


0.072 


Sample IV 


1.887 


10.9° 


< 10-3 


< 10-3 


0.072 



Table 6: For each of the four samples, the theoretical values of the amplitude and the phase 
are calculated using the response functions for (a) AGASA and (b) Haverah Park. We also 
give the probabilities that the experimental amplitude and phase come from within the 95% 
confidence limits of the model. 
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(a) Auger I 



Model 




^(^<s) 






(95%) 


(95%) 


(95%) 


Sample I 


< 10-3 


0.801 


< 10-3 


Sample II 


< 10-3 


0.811 


< 10-3 


(b) Auger II 


Model 


P{as, cr^) 




P(^^) 




(95%) 


(95%) 


(95%) 


Sample 1 


< 10-3 


0.001 


< 10-3 


Sample II 


< 10-3 


0.001 


< 10-3 



(c) Predictions for South Auger 



Model 


Amplitude 


Phase 


Amplitude 


Phase 




(hi) 


ih) 


ih2) 




Sample 1 


0.471 


49.4° 


1.187 


22.4° 


Sample 11 


0.469 


48.7° 


1.178 


22.3° 


Sample 111 


0.568 


123.2° 


1.395 


85.5° 


Sample IV 


0.411 


116.0° 


1.268 


83.9° 



Table 7: For each of the Samples, this shows the impact of the South Auger Observatory. 
Table (a) assumes the amplitude and phase of the AG AS A experiment, table (b) the ampli- 
tude and phase from Haverah Park; in both cases, 1000 events are analysed. In Tables (a) 
and (b), hues are omitted if all the probabilities vanish. Table (c) presents the predictions 
for our four Samples using the two possible response functions hi{5) and h'2,{5). 
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Fig. 8. — Plots of the amplitude against the phase for the four Samples of galaxies. Also 
shown are the 95% confidence limits for (a) the AGASA and (b) the Haverah Park exper- 
iments. The lower panels show the expected impact of South Auger. Marked are the 95% 
confidence limits assuming the amplitude and phase of (c) the AGASA and (d) the Haverah 
Park experiments, but with 1000 events. 
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4.2. Results 

To take into account the declinations of the galaxies, we use the discrete analogue-s of 
equations and (|^): 

N = J2HSi)\: (29) 

«=|em^o^, b=^j:mf^, (30) 

where the index i runs along all the galaxies in the Sample and h{6) is the response function 
of each experiment. 

Figure |^ shows the distribution of galaxies in the four Samples. The supergalactic plane 
is clearly visible, as is the Zone of Avoidance within ~ 15° of the Galactic plane. Some 
of the prominent local features are also readily apparent. The Virgo cluster is centered at 
a ~ 12^30™, 6 ~ 10°, while the Fornax cluster is at a ~ 3^, 6 ~ —40°. The densest part of the 
Local Supercluster and the Hydra-Centaurus supercluster are prominent at a ~ 12^^ — 1*^, 6 ~ 
—40°, while the Perseus-Pisces cluster stretches from a ^ 1^ — 3^,6 ^ 40°. Of course, 
any optically selected survey is strongly affected by extinction at low Galactic latitudes. 
We could avoid this bias by using a different catalogue such as the IRAS galaxies. This 
though merely substitutes one bias for another, as early-type galaxies are under-represented 
in the IRAS survey, whereas late-type dusty disk galaxies are over-represented. How serious 
is the bias caused by low latitude extinction on our calculations? Let us note that the 
Dwingeloo Obscured Galaxy Survey conducted a shallow search in the entire Northern Zone 
of Avoidance for nearby and/or massive galaxies and this yielded just 5 candidates [^. They 
estimated that there are ~ 100 galaxies missing within ~ 50 Mpc. The distance factors in 
eqs. (^,^) mean that the nearest galaxies have the largest weight. Hence, this suggests 
that low latitude obscuration will not have a deleterious effect on our results. Only if there 
are large numbers of missing galaxies within 10 Mpc is our calculation likely to be in error, 
and this circumstance does not seem to be the judged by the available evidence. 

Table ^ shows the signal expected for AGASA and Haverah Park. For all four Samples, 
the size of the amplitude is significant, so that detection of the anisotropy is not enough to 
implicate decaying dark matter haloes. Samples I and II, which contain almost all the nearby 
galaxies, do yield an anisotropic signal of the same order of magnitude as seen by AGASA 
and Haverah Park. In fact, the direction of the signal is well within the 95% confidence 
limits for the Haverah Park dataset. The probabilities that the AGASA and Haverah Park 
datasets are drawn from within the 95% confidence limits of the model are moderate (0.264) 
and high (0.858) respectively. The phase of the anisotropy is robust against changes in the 
cut-off radius, suggesting that it is the very closest galaxies (within ~ 20 Mpc) that are 
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contributing most of the signal. The direction of the signal shows that it is is dominated 
by nearby structures in the supergalactic plane, especially the Virgo cluster. Of all our 
models, Samples III and IV, which contain the nearby bright galaxies within 50 and 100 
Mpc respectively, give the largest amplitude for the anisotropy signal. They are clearly 
disfavored by the existing AGASA and Haverah Park datasets. Table |^ shows there is 
vanishingly small probability that either dataset is drawn from within the 95% confidence 
limits of the model. The observed distribution of UHECRs seems to be more isotropic than 
expected if they emanate just from the nearby bright galaxies. 

Table |^ shows the potential effects of the Auger Observatory. Again, Auger I has the 
same amplitude and phase as recorded by AGASA, while Auger II is the same as that 
recorded by Haverah Park, with the number of events increased to 1000. Of course, if 
UHECRs originate in nearby galaxies, then the Northern hemisphere signal will be different 
from the Southern hemisphere. We are simply resizing the experimental box to match 
Auger's sensitivity as an illustration of its ability to discriminate. Tables |^ (a) and (b) are 
composed mainly of zeros. The large number of events now enables all four Samples to be 
ruled out, as the signals do not coincide with those measured by AGASA and Haverah Park. 
Figure ^ shows these results pictorially, together with the 95% confidence limits. 

The signal seen by South Auger is recorded in Table |^ (c). The amplitude is greater than 
that seen by the Northern Hemisphere detectors, and the phase points in a different direction. 
For Samples I and II, we obtain results that depend on the assumed detector response so 
it is difficult to make definite predictions. The number of events for a 5a detection in the 
amplitude is ~ 120 (for the hi response) and ~ 40 (for the /i2 response). For Samples HI and 
IV, we obtain large amplitudes irrespective of the response function. Only ~ 30 events are 
needed for a 5cr detection of the amplitude. The direction of the phase is largely controlled 
by the Fornax cluster, which is the nearest rich cluster to us. Although it is less massive 
than the Virgo cluster, its effect is dominant in the Southern sky. Suppose that we wished 
to ensure that the 3a (or 5a) error bars on the measured phase at South Auger are ±20°. 
This is sufficient to implicate an extragalactic origin of the UHECRs and requires ~ 350 (or 
~ 1000) events. 



5. Conclusions 

The origin of the ultra-high energy cosmic rays (UHECRs) is unknown, but the dis- 
tribution of arrival directions provides important clues. We have examined two possible 
hypotheses for the origin in detail. As the available dataset of UHECRs is too small to 
yield definitive evidence for anisotropy, we have concentrated on the prospects for the Auger 
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Observatory, especially its southern station at Malargiie, Argentina. 

First, UHECRs may be produced by the decay of dark matter in the halo of the Galaxy. 
The offset of the Sun from the Galactic Centre causes an anisotropy signal. The magnitude 
of this anisotropy is largest for dark halo models which are cusped (such as the Navarro- 
Frenk- White or NEW profile). However, there is overwhelming astronomical evidence that 
the Galactic halo does not have a cusped density profile. It is more likely that the halo 
has a core radius ~ 10 kpc, so that dark matter is dynamically unimportant in the central 
parts but dominates more and more beyond the Solar circle. The halo may be isothermal or 
triaxial or possibly even tilted with respect to the Galactic plane. If such models describe the 
distribution of decaying dark matter, then the amplitude of the anisotropy 5 is ^ 0.4 for the 
AGASA and Haverah Park detectors, whilst the phase of the anisotropy points towards the 
direction of the Galactic Centre. The shape of the halo controls the phase, whereas the scale 
of the halo controls the amplitude. Spherical halo models yield anisotropies whose phase 
coincides with the direction of the Galactic Centre, whereas triaxial models give angular 
deviations of ~ 25° depending on the halo shape and detector location. The amplitude 
of the anisotropy is controlled by the halo core radius Tc- This parameter is not easy to 
determine by astrophysical means. So, should this mechanism of UHECR production be 
proved correct (e.g., through confirmation of the expected energy spectrum [0), then it will 
provide a direct way of measuring the core radius. 

If UHECRs originate in dark haloes, then it has been suggested that the halo of M31 
would be visible as a hotspot [T^. However, the size of this effect seems to be modest. Partly 
this is because recent analyses have revised the overall mass of the M31 halo downwards, 
making it roughly comparable to the mass of the Galactic halo [^]. Eurthermore, UHECRs 



emanating from M31 are deflected by a few degrees and so the effect must be sought in a field 
size of ~ 10° X 10°, which dilutes the signal |2^. Our calculations suggest that only a mild 
enhancement is expected (^ 60%) and its identification requires considerably more events 
than in the AGASA and Haverah Park datasets, in agreement with |^D|] but not with ^\. 



Note that the South Auger detector will not help much in this regard, as M31 is not visible 
from the latitude of Malargiie. 

A robust prediction of the decaying dark matter hypothesis is that the amplitude of 
the anisotropy at the South Auger station will be larger than at the Northern Hemisphere 
sites, simply because the Galactic Centre lies at southern declinations. Isothermal haloes, 
whether spherical, triaxial or tilted, yield anisotropies of amplitude ~ 0.5 for a halo core 
radius ~ 10 kpc. Typically, the detection of between 150 and 500 events at South Auger 
will be required for a 5a identification of the anisotropy signal. However, if the core radius 
is smaller (^ 1 kpc), then the amplitude of the anisotropy is ^ 1.0 and perhaps as few as 
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40 events would suffice. The phase will give information on triaxiality or tiltedness, and this 
may be obtained with ~ 2000 events for isothermal-likc models. More importantly, ~ 500 
events are sufficient to confirm that the phase points in the rough direction of the Galactic 
Center, which would implicate a Galactic origin of the UHECRs. 

A second hypothesis is that the UHECRs may originate in the nuclei of nearby galaxies, 
perhaps produced by supermassive black holes. If only the nearby bright galaxies are the 
sources, then the amplitude of the anisotropy is much greater than observed in the AGASA 
and Haverah Park datasets (typically 5 ~ 1.9). It therefore seems that this hypothesis can 
already be ruled out, as the probabilities that the AGASA and Haverah Park datasets are 
consistent with such signals are very small. However, some caution is needed as the number 
of UHECRs is still small and the effects of low latitude obscuration have been neglected in 
the analysis. Furthermore, the expected anisotropy may be diluted by up to a factor of ~ 2 
by a possible isotropic background. 

If extragalactic sources are to provide the UHECRs, then the population must be larger 
than just the nearby bright galaxies. The smallest galaxy presently suspected of possessing 
a supermassive black hole is M32. Samples of nearby galaxies brighter than M32 yield an 
anisotropy signal S ~ 0.5 in amplitude and ~ 120° in phase. This is in good agreement 
with the signal seen by Haverah Park, and in rough agreement with that seen by AGASA. In 
the decaying dark matter hypothesis, the phase always points in the approximate direction 
of the Galactic Centre. If, however, all nearby galaxies provide the UHECRs, then the 
phase at the AGASA and Haverah Park detectors points towards a ^ 9^ and is controlled 
by prominent mass concentrations in the supergalactic plane, such as the Virgo cluster 
(a fa 12^30"^, 5 f« 10°). 

The measurement of significant anisotropy with the South Auger station does not by 
itself validate the decaying dark matter hypothesis, as this phenomenon occurs in our ex- 
tragalactic samples as well. When all galaxies intrinsically brighter than M32 are included, 
the amplitude of the anisotropy is 5 ~ 0.6 irrespective of whether the cut-off is 50 or 100 
Mpc. The number of events for a 5a detection in the amphtude is ^ 120. The phase ■»? 
does not point toward the Galactic Centre, but is controlled by the nearby Fornax cluster 
{a fs 3**, 5 ~ —40°). Suppose that we wished to ensure that the 3a error bar on the measured 
phase at South Auger is ±20°, which would be strong evidence for an extragalactic origin. 
This requires ~ 350 events. 

In conclusion, it seems that one of the most important contributions that the South 
Auger experiment can make is to identify the phase direction in the Southern hemisphere. 
If UHECRs have a Galactic origin, then the phase will point towards the Galactic Centre. 
If UHECRs have an extragalactic origin, then it seems they must be produced by almost all 
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nearby galaxies (else the signal seen by AG ASA and Haverali Park would be much larger). 
Then, the phase of the anisotropy recorded by South Auger will lie roughly in the direction of 
the Fornax cluster. Provided the effects of the intergalactic magnetic field can be neglected, 
this gives an unambiguous discriminant between the two theories and requires only ~ 350 — 
500 events. This should be obtained within the first three years of operation of South Auger. 
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